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Abstract 

A hydrodynamic model coupled to the statistical hadronization code Therminator is used to study 
a set of observables in the soft sector at RHIC. A satisfactory description of the /? ± -spectra and 
elliptic flow is obtained, similarly to other hydrodynamic models. With the Gaussian initial 
conditions the transverse femtoscopic radii are also reproduced, providing a possible solution of 
the RHIC HBT puzzle. 



1. Intoduction 

The consistent description of various features of the soft hadron production in the nucleus- 
nucleus collisions at the Relativistic Heavy-Ion Collider (RHIC) is a well known problem [1]. 
The so called RHIC HBT puzzle |fl I! H II refers to the difficulty of simultaneous description 
of the hadronic transverse-momentum spectra, the elliptic flow coefficient V2, and the Hanbury- 
Brown-Twiss (HBT) interferometry data in various approaches including hydrodynamics J31 la 

bus. 

Recently, we have constructed a hydrodynamic code ifloilTlll that is coupled to the statistical 
hadronization model Therminator [12]. Within this approach a successful uniform description of 
the soft hadronic RHIC data has been accomplished 111 311 . The main ingredients of our approach 
are the following: i) We use a realistic equation of state that interpolates between the lattice 
QCD results at high temperatures and the hadron-gas results at lower temperatures CI. ii) The 
single freeze-out scenario including all well established resonance states is assumed 11511 . This 
assumption leads to shorter emission times and helps to reproduce the ratio R out /Rside- iii) The 
use of Therminator allows for the use of two-particle methods in the evaluation of the correlation 
functions (with and without the Coulomb corrections), iv) The early starting/thermalization time 
Tj = 0.25 fm is assumed. This helps to develop fast the transverse flow and shortens the evolution 
time, v) The initial conditions for the energy density in the transverse plane are taken in the 
Gaussian form lfl3ll 

I x 2 y 2 \ 

e (xJ = £l exp^— -— j. (1) 
The values of the width parameters a and b depend on the centrality class The Gaussian profiles 



are obtained with the Monte-Carlo Glauber simulations done with Glissando II 1611 . vi) The initial 
conditions include fluctuations of the initial eccentricity I^, lr7ll . 
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2. Results 



In Fig.Q]we show our model results for the transverse-momentum spectra and compare them 
with the PHENIX data HH for Au+Au collisions at the top RHIC energy of -y/^wv = 200 GeV. 
One observes a very good agreement between the model predictions and the data. The model 
results have been obtained with the initial central temperature Tj = 460 Me V and the final/freeze- 
out temperature Tf = 145 MeV. These are the two main parameters of our approach. The values 
of the width parameters follow from the Monte-Carlo Glauber simulations and for the centrality 
class 20-30% we have found a = 2.00 fm and b = 2.59 fm. 




Figure 1: The transverse-momentum spectra of pions, kaons and protons for the centrality class 20-30%. The model 
results with the Gaussian initial conditions are plotted as functions of the transverse momentum and compared to the data 
from O. 
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Figure 2: The elliptic flow coefficient V2 for the centrality class 20-40%. The model results with the Gaussian initial 
conditions are plotted as functions of the transverse momentum and compared to the data from [19]. Model parameters 
are the same as in Fig.[T] 

In Fig.|2]we show our model results for the elliptic flow coefficient V2. The same values of the 
parameters have been used as in the calculation of the spectra. The pion+kaon data are very well 
reproduced, while the model results for the protons slightly overshoot the data. This behavior 
may be attributed to the lack of hadron rescattering in the final state in our approach. 

Our HBT results are shown in Fig. [3] Again, the same values of the input parameters have 
been used. We observe good agreement between the data [20] and the model calculations. In 
particular, the ratio R ou t/Rside is very well reproduced. 
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Figure 3: The pion HBT radii R^ e , R mt , R\ om , and the ratio /f ut/^side f° r the centrality class 20-30%. The 
results of the model calculation with the Gaussian initial conditions (lines) are compared to the data from 
J2(ill (dots). Model parameters are the same as in Fig.[TJ 

Encouraged by the success of reproducing the azimuthally averaged HBT radii, we have 
also calculated the azimuthally sensitive femtoscopic observables for different centralities and 
average transverse momenta 112 ill . The summary of our results is shown in Fig. [4] The model 
results are compared to the experimental STAR data J22ll . For each centrality, associated here 
with the number of participants N p!irt on the horizontal axis, we plot the experimental points 
(filled dots) and the model results (empty symbols). The points from top to bottom correspond 
to k T contained in the bins of 0.15-0.25 GeV, 0.25-0.35 GeV, and 0.35-0.6 GeV. The top panels 
show the radii squared averaged over the <p angle, from left to right, R^ ulQ , ^jeQ' m ^ ^longer 
The bottom panels show the magnitude of the allowed oscillations divided by /ft . () , which is 
the adopted convention used in presenting the experimental data. 

We conclude with the statement that the consistent and uniform description of the soft hadronic 
data at RHIC may be achieved within the hydrodynamic approach if a proper choice of the initial 
profile and a realistic equation of state are used. We note that a similar conclusion has been also 
reached recently in Ref. Ir23ll (see also J3) where the fast building of the transverse flow that is 
required for the correct description of the HBT radii is achieved with the very early start of the 
hydrodynamics (t; =0.1 fm) and the inclusion of the viscous effects. However, the V2 coefficient 
is not evaluated in lr23ll . 
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Figure 4: Results for the RHIC HBT radii and their azimuthal oscillations. For each value of Np M on the 
horizontal axis we plot the experimental points (filled symbols) and the model results (empty symbols). The 
points from top to bottom at each plot correspond to k T contained in the bins 0. 15-0.25 GeV (circles), 0.25- 
0.35 GeV (squares), and 0.35-0.6 GeV(triangles). The top panels show R 2 oM0 , ^sideO' an ^^?ong0' tne bottom 



panels the magnitude of the allowed oscillations divided conventionally by R 
Model parameters for different centralities are given in (2lll . 
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